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I.  INTRODUCTION 

This  memorandum  describes  two  versions  of  a  program  for  analyzing  linear, 
passive  electrical  circuits.  The  first  version  is  analytic  and  calculates 
steady  state  frequency  response  while  the  second  is  numerical  and  calculates 
complete  response  to  initial  conditions  and/or  forcing  functions.  Both 
versions  prompt  the  user  to  describe  the  circuit  by  specifying  the  number  of 
nodes  and  the  components  which  interconnect  them.  Any  interconnection  of 
resistors,  capacitors  and  inductors  is  allowed.  The  steady  state  program,  and 
the  complete  response  program  with  a  forcing  function,  are  configured  for  two 
port  circuits.  A  common  ground  is  assumed. 

The  steady  state  version  requires  that  the  frequency  range  of  interest  be 
specified.  In  addition,  the  admittance  of  the  driving  source  must  be 
specified.  In  the  complete  response  version  the  user  must  choose  a  time  step 
appropriate  to  the  problem  being  considered.  The  u3er  is  prompted  to  supply 
the  necessary  initial  conditions  and  is  offered  a  selection  of  forcing 
functions. 

The  programs  are  written  in  basic.  The  program  list  for  the  steady  state 
version  is  given  in  Appendix  1  and  that  for  the  complete  response  version  in 
Appendix  2.  The  programs  are  available,  for  the  IBM  PC,  on  5  %  "  floppy  disk. 

II.  THE  STEADY  STATE  FREQUENCY  RESPONSE  PROGRAM 

To  use  the  program,  the  circuit  nodes  are  numbered  starting  with  number 
one  at  the  input  and  proceeding  to  the  highest  number  at  the  output.  The 
ordering  for  the  intermediate  nodes  is  arbitrary  and  the  ground  node  is  not 
numbered.  The  program  assumes  that  a  sinusoidal  current  source  is  attached 
between  node  one  and  ground.  The  program  sets  the  strength  of  this  current 
source  so  as  to  combine  with  the  specified  shunt  conductance  to  yield  a  unit 
voltage  source  (open  circuit).  The  values  of  the  circuit  components  are 
entered  as  admittances  (i.e.,  as  capacitance  in  farads,  conductance  in  mhos, 
and  inverse  inductance  in  inverse  Henrys).  By  this  means  the  absence  of  any 
particular  component  is  indicated  by  inputing  the  value  zero.  For  example  the 
absence  of  a  resistor  connecting  two  nodes  is  indicated  by  specifying  zero 
conductance. 
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A  lowest  frequency  (FI  in  the  program,  here  called  f)  is  entered  (in  Hz 
and  as  a  whole  power  of  ten)  together  with  the  number  of  decades  to  be 
covered.  The  output  node  voltage  and  its  phase  relative  to  the  input  are  then 
tabulated  for  the  frequencies  f,  2f,  5f.  10f,  20f,  50f,  . 

An  example  of  the  use  of  the  steady  state  program  is  given  in  Appendix 
3.  Some  useful  equivalencies  are  given  in  Appendix  5. 

III.  THEORY  FOR  THE  STEADY  STATE  PROGRAM 

In  this  program  the  variables  are  the  node  voltages.  Conductance  is 
symbolized  by  G,  capacitance  by  C,  and  the  inverse  of  inductance  by  B. 
Subscripts  on  these  parameters  indicate  the  nodes  which  the  component 
interconnects.  Thus  Gy  represents  a  conductance  between  node  i  and  node  j. 

A  repeated  index  is  used  to  represent  a  connection  between  a  node  and 
ground.  Thus  G^  is  the  conductance  between  node  i  and  ground.  If  we  sum  the 
currents  into  each  node,  assume  all  variables  to  have  a  eia)t  time  dependancy, 
separate  each  voltage  into  its  real  and  imaginary  parts,  and  rearrange  terms, 
there  results  a  set  of  simultaneous  equations  as  given  by  Equation  1 .  A 
current  source  of  zero  phase  and  strength  I  is  assumed  to  be  attached  between 
node  1  and  ground.  In  the  program,  lines  300  through  350  incorporate  the 
"self"  terms  into  the  diagonal  sub-matrices.  Lines  490  through  720  accumulate 
the  sums  in  the  diagonal  sub-matrices  and,  in  addition,  fill  in  the  off 
diagonal  sub-matrices.  The  restriction  of  the  index  J  to  values  greater  than 
I  (line  370)  is  an  economy  permitted  by  the  identity  of  elements  with 
interchanged  indices,  i.e.,  G^2  *  &2\  e^c* 

The  frequencies  for  which  these  equations  are  to  be  solved  are 
established  by  lines  880  through  960.  In  the  subroutine  represented,  by  lines 
980  through  1380  the  C,  G,  and  B  matrices  are  summed  for  each  frequency  (line 
1010)  the  resulting  simultaneous  equations  solved  using  a  Gauss-Jordan 
elimination  method  (lines  1080  through  1280)  and  the  phase  and  amplitude  at 
the  output  node  calculated  (lines  1290  through  1350). 


IV.  THE  COMPLETE  RESPONSE  PROGRAM 

The  circuit  description  is  entered  in  this  program  in  the  same  manner  as 
in  the  steady  state  program.  That  is,  by  specifying  the  number  of  nodes  and 
the  components  interconnecting  them.  When  a  forcing  function  is  specified,  it 
is  assumed  to  be  attached  between  node  one  and  ground.  The  output  is  the 
highest  number  node.  A  time  step  and  number  of  steps  appropriate  to  the 
problem  under  consideration  must  be  entered.  The  time  step  must  be  small 
compared  to  the  time  associated  with  any  oscillation  or  decay  within  the 
circuit. 

An  example  of  the  use  of  the  complete  response  program  is  given  in 
Appendix  4. 


V.  THEORY  FOR  THE  COMPLETE  RESPONSE  PROGRAM 

By  summing  the  currents  into  each  node  and  re-arranging  terms  we  obtain 
an  equation  for  each  node  of  the  form. 


Z  +  V(I.a)  £  G(I,J)  *  (  /Vl.aOdt)  Z  B (I,J) 


-  Z  C(l,j)  -  Z  V(J,a)  G( I, J)  -  E  (/  V( J.aOdt]  B(I,J)  ♦  IIO(I) 

J*I  Jtl  J*I  o 


where : 


V( I, a)  -  voltage  at  node  I  for  either  the  present  or  the  next  step 

C(I,J)  -  capacitance  connecting  node  I  to  node  J 

G(I,J)  -  conductance  connecting  node  I  to  node  J 

B(I,J)  *  inverse  inductance  connecting  node  I  to  node  J 

1(1)  -  current  sources  attached  between  node  I  and  ground 
o  n  o 

IIO(I)  -  (J  VU.cOdt)  £  B( I, J)  -  £  (J  V(J,a)dt)  B(I, J) 

J-1  Jtl  — 


V. V. 


and  "self"  terms  (e.g.,  C(I,I))  represent  connections  to  ground.  In  the 
program  the  following  additional  variables  are  employed: 


D(I,a) 

S(I,a) 

C(I)  - 

G(I)  - 

B(  I)  - 

DS( I)  - 

VS(I)  - 

SS(I)  - 


avd.tt) 

3t 


-  J  V(I,a)dt 

o 

N 

I  C(I,J) 
J-1 

N 

Z  G( I , J) 

J-1 

n 

Z  B(I, J) 
J-1 


Z  D( J,a)  C(I,J) 
J-I 


£  V(J,a)  G(I, J) 
J*I 


Z  S(J,a)  B( I , J) 
J*I 


A  second  index,  a  -  1  or  2,  has  been  added  to  the  voltage  to  indicate 
either  the  present  time  step  or  the  next  time  step  forwards.  This  has  been 
done  so  that  the  equations  can  be  integrated  using  a  variation  of  a  second 
order  Runge-Kutta  method. 

Lines  320  through  860  assemble  the  description  of  the  circuit  (i.e., 
evaluate  C(I,J),  G(I,J)  and  B(I,J)  for  all  I  and  J)  and  establish  the  initial 
conditions  (evaluate  V(I,1)  and  IIO(I)). 

The  introduction  of  forcing  functions  is  accomplished  through  lines  870 
through  1110.  When  an  abrupt  change  is  made  in  the  current  1(1)  entering  a 
node,  several  steps  of  simple  (Eulerian)  integration  will  re-establish  self- 
consistant  values  of  the  various  derivatives.  The  second  order  Runge-Kutta 
does  not  absorb  large  inconsistancies  and  progress  stably.  Therefore  the 
calculation  starts  with  several  steps  of  Eulerian  integration  using  a  very 
small  time  step  (subroutine  1880). 
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The  use  of  Equation  2  to  construct  an  algorithm  is  most  clearly  seen  in 
subroutine  1880.  For  any  given  node  not  all  the  terms  in  Equation  2  may  be 
present.  This  occurs  when  some  types  of  components  are  absent,  (e.g.,  for  a 
node  having  no  capacitve  connection,  (C(I)  -  0).  Lines  1980  and  2030  select 
the  appropriate  form  of  Equation  2  to  give  the  highest  voltage  derivative 
necessary  for  each  node.  It  is  then  possible  to  perform  the  Eulerian 
integration  to  move  forwards  in  time.  Note  that  for  a  node  joining  only 
inductors,  that  no  integration  is  necessary  because  only  S(I,a)  needs  to  be 
known  and  this  is  found  algebraically  from  Equation  2  (line  2070).  Similarly 
only  a  single  integration  (line  2050)  is  necessary  for  a  node  without 
capacitors  because  only  V(I,a)  and  S(I,at)  need  to  be  known. 

The  main  algorithm  is  contained  in  subroutine  1220.  An  initial 
integration  is  performed  in  lines  1230  through  1300.  A  first  estimate  of  the 
necessary  derivatives  at  the  next  time  step  is  made  in  lines  1310  through 
1480.  This  estimate  is  made  using  those  updated  values  which  the  initial 
integration  made  available.  The  integration  is  then  repeated  in  lines  1480 
through  1550  using  the  average  of  tne  derivatives  at  the  present  and  next  time 
steps.  In  lines  1560  through  1720  the  estimate  of  the  derivatives  at  the  next 
time  step  is  refined  by  using  the  improved,  updated  values.  A  last  repetition 
of  the  integration  is  performed  in  lines  1730  through  1800. 


. 
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VI.  CONCLUSION 

The  author  finds  great  utility  in  having  a  program  that  establishes  and 
simplifies  the  simultaneous  equations  for  arbitrary  circuit  configurations. 

To  write  these  equations  by  hand,  for  each  specific  case,  has  always  been  a 
tedious  and  error  prone  process. 

Among  the  problems  to  which  the  program  has  been  applied  are,  the  design 
of  the  fast  passive  integrator  (NRL  MR  4939),  the  design  of  a  pulse  forming 
capacitor  bank,  a  lumped  parameter  modeling  of  the  field  problem  for  an  e-beam 
in  an  air  filled  drift  tube,  the  analysis  of  several  electro-accoustic  systems 
and  the  verification  of  several  5-pole  Butterworth  filter  designs.  All  of 
these  problems  involve  circuits  of  S  10  nodes  for  which  an  analysis  by  hand 
would  be  prohibitively  long  yet  the  analysis  using  a  desk  top  computer 
requires  only  a  few  minutes. 
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APPENDIX  1 


10  REM  Steady  State  Program 

20  REM  This  program  is  con-figured  far  a  two  port  circuit.  Nodes  are  numbered 
from  the  input  which  is  number  one,  to  the  output  which  is  number  N. 

30  REM  The  ground  is  not  numbered.  A  Source  of  one  volt  (open  circuit)  and  of 
the  specified  admittance  is  assumed  to  be  attached  to  node  one. 

40  REM  Resistance  is  in  Ohms;  capacitance  in  Farads;  inductance  in  Henries; 

frequency  in  Herts;  and  the  output  amplitude  is  in  Volts. 

50  LPRINT  TAB (5) "Steady  state  circuit  analysis  program,  version  1“ 

60  LPRINT  TAB (5)  DATE* 

70  LPRINT 

30  INPUT  "no.  of  nodes  7";N 

90  N1=2*N 

100  N2=N1+1 

110  N3=N1— 1 

120  N4=N 1 +2 

130  LPRINT  TAB (5) "Steady  state  analysis  for  a  circuit  with  ";N;  "nodes" 

140  LPRINT 

150  DIM  G(N4.N4) 

160  DIM  C(N4,N4) 

170  DIM  B(N4,N4) 

180  DIM  A (N4, N4) 

190  FOR  1=1  TO  N 

200  PRINT  "1/R,C, 1/L  7  node  " ; I ;  "to  gnd” 

210  INPUT  G 
220  INPUT  C 
230  INPUT  B 

240  LPRINT  TAB(5)"Node  " ;  I ;  "  has  1/R  to  ground  of  ";G 

250  LPRINT  TAB (3) "  C  to  ground  of  ";C 

260  LPRINT  TAB (3)"  1/L  to  ground  of  ";B 

270  LPRINT 
230  L=2*I-1 
290  P=L+1 

300  C(L,L)=C(L,L)+C 
310  B<L,L)=B(L,L>— B 
320  C<P.P)=C(P,P)+C 
330  B  <P, F) -B (P, P) -B 
340  G(L.P)=G(L,P)+G 
350  G(P,L)=G(P.L)-G 
360  IF  I=N  THEN  GOTO  750 
370  FOR  J=I+1  TO  N 
330  M=2*J-1 
390  Q=M+1 

400  PRINT  "l/R.C.l/L  7  node  ";I;"to  node  ";J 
410  INPUT  G 
420  INPUT  C 
430  INPUT  B 

440  LPRINT  TAB(5)"From  node  ”;I;”  to  node  ";J 


450 

LPRINT 

TAB (5) " 

1/R 

460 

LPRINT 

TAB (5) " 

C 

470 

LPRINT 

TAB (5) " 

1/L 

430 

LPRINT 

490 

C  (L,M>  • 

=-C 

500 

B(L,M) 

=B 

510  C (P, Q) =— C  . 

520  B  <F',  Q)  =B  \ 

530  C  (M,  Li  =— C  '! 

540  B(M,L)=B 
550  C(Q,F)=-C 
560  B (Q, P ) =B 
570  GCL.Q)=-G 
580  G(P,M)=G 
590  S(Q,L)=B 
600  G(M,P)=-G 
610  C(L,L>=C(L,L)+C 
620  B  <L, L)  =B  <L,  L)  -B 
630  C(M,M)=C(M,M)+C 
640  B<M,M)=B(M,M)-B 
650  C(P,Pi=C<P,P)+C 
660  B(P,P)=B(P,P)-B 
670  C(Q,Q)=C(D,Q)+C 
6 SO  B  (Q,  Q)  =B  <Q,  Q)  — B 
690  G(L,P)=G(L.P)+G 
700  G(P,L)=G(P,L>-G 
710  G(M,D)=G(M,Q)+G 
720  G(Q,M)=G(Q,M)-G 
730  NEXT  J 
740  NEXT  I 

750  PRINT  “Source  admittance  ?” 

760  INPUT  IC 

770  LPRINT  TAB (5) "Source  admittance  =  ";IC 
780  G  <  1 , 2 ) =G ( 1 , 2 ) + 1 C 
790  G  <2, 1 ) =G (2, 1 ) - 1 C 

800  PRINT  "Low  -frequency,  no.  o-f  decades  ?" 

810  INPUT  FI 
820  INPUT  D 

830  LPRINT  USING  "  Lowest  -frequency  considered  #.### 

840  LPRINT  TAB(5)"No.  of  decades  considered  ";D 
850  LPRINT 

860  PRINT  TAB  < 1 ) "Freq. ";TAB(12i "Phase" ; TAB (25) "Ampl . " 

870  LPRINT  TAB (5) "Freq. ”;TAB(16) "Phase" ; TAB ( 25 ) "Ampl  .  " 

880  FOR  K=0  TO  D-l 
890  F=F 1  *  10 'K 
900  FOR  L= 1  TO  2 
910  W=6. 283185*L*F 
920  GOSUB  980 
930  NEXT  L 
940  W=31 . 415927#*F 
950  GOSUB  980 
960  NEXT  K 
970  GOTO  1390 
930  FOR  1=1  TO  Ml 
990  FOR  J=1  TO  N1 
1 000  A ( I , J  >  =0 

1010  A ( I , J) =W*W*C  < ;  I , J) +W*G ( I ,  J) +B ( I , J) 

1020  NEXT  J 
1030  NEXT  I 


1040  FOR  1=1  TO  N1 


1050  A<I,N2  =0 

1060  NEXT  I 

1070  A<2,N2)=-IC*W 

1080  FOR  J=1  TO  N1 

1090  FOR  I=J  TO  N1 

1100  IF  A ( I , J ) <  >0  THEN  GOTO  1120 

1110  NEXT  I 

1120  FOR  M=1  TO  N2 

1130  X=A<J,M) 

1140  A  ( J ,  M) =A ( I , M) 

1150  A  < I , M) =X 
1160  NEXT  M 
1170  Y=1/A(J,J) 

1180  FOR  M-l  TO  N2 
1190  A<J,M)=Y*A<J,M) 

1200  NEXT  M 

1210  FOR  1=1  TO  N1 

1220  IF  I=J  THEN  GOTO  1270 

1230  Y=-A ( I , J ) 

1240  FOR  M=1  TO  N2 

1250  A<I,M)=A<I.M)+Y*A<J,M> 

1260  NEXT  M 
1270  NEXT  I 
1280  NEXT  J 

1290  V=SQR(A(N3,N2) *A <N3, N2) +A <N1 , N2) *A<N1, N2) ) 

1300  P=57.3*ATN(A(N1,N2) /A(N3,N2) ) 

1310  IF  A (N3, N2) >0  THEN  GOTO  1360 
1320  IF  A  <N1 , N2) >0  THEN  GOTO  1350 
1330  P=P-1S0 
1340  GOTO  1360 
1350  P=P+180 

1360  PRINT  USING  +###  W/6. 283185. INT <P+. 5) , V 

1370  LPRINT  USING  "  +###  - W/6. 283185.  INT  CP+.  5)  ,  V 

1380  RETURN 
1390  END 
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at  nodes" 


wi th  "  ;  N; 


10  REM  Complete  Response  Program 

20  REM  Nodes  are  numbered  -from  the  input  which  is  number  to  the  outpu 
number  N.  The  ground  is  not  numbered. 

30  REM  Resistance  is  in  Ohms:  capacitance  in  Farads;  inductance  in  Henries;  ti 
in  seconds;  and  the  output  in  Volts. 

40  LPRINT  TAB<5)  "Complete  responce  circuit  analysis  program,  version 
50  LPRINT  TAB (5)  DATE* 

60  LPRINT 

70  PRINT  "No.  at  nodes" 

80  INPUT  N 

90  LPRINT  TAB (5) "Compl ete  responce  analysis  tor  a  circuit  with  ";N;  " 

100  LPRINT 

110  PRINT  "Time  step" 

120  INPUT  TS 

130  PRINT  "No.  of  steps" 

140  INFUT  M 

150  LPRINT  TAB<5)"time  step  =  “;TS;  "  ,  number  ot  steps  -  ";M 

160  LPRINT 
170  DIM  I (N) 

180  DIM  D (N, 2) 

190  DIM  V ( N , 2 ) 

200  DIM  S (N, 2) 

210  DIM  1 0  < N , N ) 

220  DIM  C(N,N) 

230  DIM  G  ( N ,  N ) 

240  DIM  B(N,N) 

250  DIM  CC  < N ) 

260  DIM  GG(N) 

270  DIM  BB  < N ) 

2 SO  DIM  DS  <N) 

290  DIM  VS  <M> 

300  DIM  SS(N) 

310  DIM  II0(N) 

320  FOR  1=1  TO  N 

330  PRINT  "1/R,C,1/L  tram  node  ";I;"  to  Gnd. " 

340  INPUT  G 
350  INPUT  C 
360  INPUT  B 

370  LPRINT  TAB (5) "node  ";I;  "has  1/R  to  ground  ot  ";G 
380  LPF'INT  TAB<5)"  C  to  ground  ot  ";C 


nodes. " 


390  LPRINT  TAB (5)" 

400  LPRINT 
410  G  < I . I  1 =G 
420  C  I ,  I  )  =C 
430  B ( I , I ) =B 

440  IF  B=0  THEN  GOTO  490 
450  PRINT"  Ini  ti  a),  current 
460  INPUT  10(1,1) 

470  LPR I NT  TAB (5 ) " I n l t l al 
10(1, I>;  "  amp s . " 

480  LPRINT 

470  IF  I =N  THEN  GOTO  740 
500  FOR  J =1+1  TO  N 


C  to  ground  ot  ";C 
1 /L  to  ground  ot  " ; B 


into  node"; I; "via  the  inductor  to  gnd." 


j  nto  node 


"via  the  inductor 


gnd.  is 


'  to  node  ";J 
“;G 


720  NEXT 
730  NEXT 


510  FRINT  "l/R.C.l/L  from  node  ";I;"  to  node  ";J 
520  INPUT  G 
530  INPUT  C 
540  INPUT  B 

550  LPRINT  TAB (5) "From  node  " ;  I ;  "  to  node  ” ; J 

560  LPRINT  TAB (5)"  1/R  =  ";G 

570  LPRINT  TAB (5) "  C  =  “ ;  C 

580  LPRINT  TAB (5) "  I/L  =  ";B 

590  LPRINT 

600  IF  B=0  THEN  GOTO  650 

610  PRINT" Initial  current  into  node"; I; "via  the  inductor  to  node";J 
620  INPUT  I0(I,J> 

630  LPRINT  TAB (5) " Ini tial  current  into  node  ” ; I ;  "  via  the  inductor  to  node 
"  is  "; 10(1,  J)  ;  "  amps. " 

640  LPRINT 

650  G ( I , J ) =G 

660  G ( J, I )  =G 

670  C ( I , J ) =C 

680  C ( J, I ) =C 

690  B ( I , J ) =B 

700  B  < J, I ) =B 

710  10 (J, I ) =-I0 ( I , J ) 

720  NEXT  J 

730  NEXT  I 

740  FOR  1=1  TO  N 

750  FOR  J=1  TO  N 

760  GG ( I ) =GG ( I ) +G ( I «  J ) 

770  CC ( I ) =CC ( I ) +C ( I , J) 

780  BB ( I ) =BB ( I ) +B  < I , J ) 

790  II0(I)=II0(I)+I0(I,J) 

800  NEXT  J 

810  IF  CC(I)=0  THEN  GOTO  860 

820  PRINT  "Initial  voltage  on  node"; I 

830  INPUT  V(I,  1) 

840  LPRINT  TAB (5) " Ini ti al  voltage  on  node  " ; I ;  "  is  ”;V(I,1);  "  volts." 

850  LPRINT 
860  NEXT  I 

870  PRINT" Input  a  choice  of  farcing  function" 

880  PRINT"0  =  no  forcing  function" 

890  PRINT" 1  =  unit  impulse" 

900  PR I NT "2  =  unit  step" 

910  INPUT  F 

920  LPRINT  TAB(5)"The  farcing  function  is  type  ";F 

930  LPRINT 

940  T=TS/ 10000 

950  IF  F=0  THEN  GOTO  1110 

960  PRINT"Source  admittance" 

970  INPUT  A 

980  LPRINT  TAB (5) "Source  admittance  =  ";A 

990  LPRINT 

1 000  G ( 1 , 1 ) =G  < 1 , 1 ) +A 

1 0 1 0  GG ( 1 ) =GG ( 1 ) + A 

1020  IF  F=1  THEN  GOTO  1060 

1030  1(1) =A 


" ; V ( I , 1 ) ;  "  volts. " 


»  •  •  *  *  »  “  %  *  ■  •  «  '  •  •*•**»■»•,* 


C 


1040  (30SUB  1890 
1050  SOTO  1120 
1 060  I < 1 > - 1 0004A/TS 
1070  GOSUB  1880 
1080  I(1)=0 
1090  GOSUB  1880 
1100  GOTO  1120 
1110  GOSUB  1380 
1120  T=T3 
1130  V  <N,  2) =V <N, 1 ) 

1140  LPRINT  TAB (5) "Ti me  Voltage" 

1150  FOR  L=0  TO  I NT  <M/10> 

1160  IF  CC ( N ) =0  AND  GG(N)=0  THEN  V <N, 2) = (S <N, 2) -S <N, 1 ) ) /T 
1170  PRINT  10*L*T,V(N,2) 

1130  LPRINT  USING  "  » ;  10*L*T,  V  <N,  2) 

1190  GOSUB  1220 
1200  NEXT  L 
1210  END 

1220  FOR  K=0  TO  9 

1230  FOR  1=1  TO  N 

1240  IF  CC(I)=0  THEN  GOTO  1280 

1250  V  < 1 , 2) =V (1,1) +D  <  I , 1 ) *T 

1260  S (1 , 2) =S ( 1 , 1 ) +V  < 1 , 1 )  *T 

1270  GOTO  1300 

1280  IF  GG  < I ) =0  THEN  GOTO  1300 

1290  S (1 , 2) =S ( 1 , 1 ) +V ( 1 ,  1 ) *T 

1300  NEXT  I 

1310  FOR  1=1  TO  N 

1 320  DS ( I ) =-D  < 1 , 1 ) *C ( I . I ) 

1330  VS  ( I ) =-V ( 1 , 1 >  *G ( I , I > 

1340  SS  < I ) =-S < 1 , 1 ) *8  < I . I ) 

1350  FOR  J=1  TO  N 

1360  DS(I)=DS(I)+D(J,  1)  *C  ( I ,  J) 

1370  vsd)=vs<i)+v<j,  i)*ea,j) 

1 380  SS (I ) =SS  < I ) +S ( J, 1 ) *B ( I , J) 

1390  NEXT  J 

1400  IF  CC<I)=0  THEN  GOTO  1430 

1410  D< 1,2)  =  (DS(Ii  +VS  < I ) +SS  ( I )  +1 (I) +110(1) -V<I,2) *GG(  I ) -S ( 1,2) *BB ( I) ) /CC( I) 
1420  GOTO  1470 

1430  IF  GG ( I ) =0  THEN  GOTO  1460 

1440  V(I,2)=(VS(I)+SS(I)+I (I)+II0(I)-S(I,2)*BB(I> )/GG(I) 

1450  GOTO  1470 

1460  S< I,2)=(SS(I)+I (I)+II0(I) ) /BB ( I ) 

1470  NEXT  I 

1480  FOR  1=1  TO  N 

1490  IF  CC(I)=0  THEN  GOTO  1530 

1500  S ( I , 2) =S (1,1) +T* (V( 1,2) +V( I , 1 ) ) /2 

1510  V( I .2) =V ( I, 1) +T* (D( I, 2) +D( I , 1 ) ) /2 

1520  GOTO  1550 

1530  IF  GG  < I ) =0  THEN  GOTO  1550 

1540  S ( I , 2) =S ( I , 1 ) +T* (V ( I , 2) +V( I , 1 ) ) /2 

1550  NEXT  I 

1360  FOR  1=1  TO  N 

1570  DS ( I ) =-D( 1,2) *C(  I , I) 


1 5 SO  VS ( I ) =— V ( I , 2 ) *G ( I , I ) 

1590  SSCD=-SCI,2>*BCI,  I) 

1600  FOR  J=1  TO  N 

1610  DS ( I ) =DS  ( I)  +D  C J , 2) *C ( I , J ) 

1620  VS ( I ) =VS <  I)  +V ( J , 2 ) *G(I, J) 

1630  SS ( I ) =SS  < I ) +S (J,2) *B  < I , J ) 

1640  NEXT  J 

1650  IF  CC(I)=0  THEN  GOTO  1680 

1660  D  C 1 , 2)  =  CDS  C I ) +VSC  I ) +SS  C I ) +1 < I ) +1  IOC  I ) -VC  1,2) *GG ( I > -3 C I , 2) *BB  CD) /CC ( I ) 
1670  GOTO  1720 

1630  IF  GG  < I ) =0  THEN  GOTO  1710 

1690  V  <  1 , 2 )  =  ( VS  C I )  +SS  ( I )  + 1  CD+IIOCI)  -S  C 1 , 2)  *8B  C I )  )/GG(I) 

1700  GOTO  1720 

1710  S(I,2)  =  CSS(D+I  CD+IIOCI)  )  /BBCI) 

1720  NEXT  I 

1730  FOR  1=1  TO  N 

1740  IF  CC  C I) =0  THEN  GOTO  1780 

1750  SCI,2)=SCI, 1 ) +T#  C V  C I , 2)  +V  ( 1 . 1) )/2 

1760  V  C 1 , 2 ) =V  (1,1) +T*  CD  C I, 2) +D ( I, 1 ) ) /2 

1770  GOTO  1800 

1780  IF  GG  C I) =0  THEN  GOTO  1800 

1790  S  ( 1 , 2)  =S  C 1 ,  1 )  +T *  C  V  <  1 , 2 )  -*-V  <  1 ,  1 )  )  /2 

1800  NEXT  I 

1810  FOR  1=1  TO  N 

1820  DC  I, 1)=DCI,2) 

1330  VC  I, 1)=V(I,2) 

1840  SCI, 1)=SCI,2) 

1850  NEXT  I 

1860  NEXT  K 

1870  RETURN 

1880  FOR  K=0  TO  9 

1890  FOR  1=1  TO  N 

1900  DS  C  D  =-D  C 1 , 1 ) *C  C I , I) 

1910  VS C I) =-V< I, 1) *GCI, I) 

1920  SS ( I) =-S (1,1 ) *BCI, I) 

1930  FOR  J=1  TO  N 

1940  DS  < I ) =DS  C I ) +D  C J , 1 ) *C ( I , J ) 

1950  VS  < I ) =VS  < I ) +V (J, 1 ) *3  C I , J) 

1960  SS ( I ) =S5  < I) +S ( J , 1 ) *B  < I , J ) 

1970  NEXT  J 

1980  IF  CC  ( I )  =0  THEN  GOTO  2030 

1990  DC  I ,  1  >  =  (DS<  I)  +VSC I )  +SSCD+I  <  I)  +1 10  C I )  -V  C 1 ,  1 )  *GG(D-S(I,  1)  *BB(  I )  )  /CCCI) 
2000  V ( 1 , 1 ) =V  C 1 , 1 ) +D  C 1 , 1 ) *T 
2010  S  <  1 ,  1 )  =S ( 1 , 1 ) +V  < 1 ,  1 )  *T 
2020  GOTO  2080 

2030  IF  GGCD-0  THEN  GOTO  2070 

2040  VC  I, 1)  =  (VS( I) +SSCD  +1 ( I) +1 10 C I) -S C 1 , 1 ) *BB< I ) ) /GGCI) 

2050  S  C 1 , 1 ) =S  C 1 , 1 ) +V ( 1 , 1 ) *T 
2060  GOTO  2080 

2070  S  < 1 ,  1)  =  CSSC I) +1 C I) +1  IOC  I) ) /BB  < I ) 

2080  NEXT  I 
2090  NEXT  K 
2100  RETURN 
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APPENDIX  3 

This  appendix  presents  the  results  of  a  steady  state  analysis  of  the 
following  circuit.  All  quantities  are  given  as  admittances.  The  node  numbers 
are  indicated  and  the  source  is  also  shown. 


Steady  state  circuit  analysis  program,  version  1 

02-01-1985 

Steady  state  analysis  for  a  circuit  with  6  nodes 

Node  1  has  1/R  to  ground  of  0 

C  to  ground  of  0 

1/L  to  ground  of  0 

From  node  1  to  node  2 
1/R  =  20 

C  =  0 

1/L  =  0 

From  node  1  to  node  3 

1/R  =  0 

C  =0 
1/L  =  0 

From  node  1  to  node  4 

1/R  =  0 

C  =  0 

1/L  =  0 

From  node  1  to  node  5 

1/R  =  0 

C  =  0 

1/L  =  0 

From  node  1  to  node  6 

1/R  =  0 

C  -  0 

1/L  =  0 

Node  2  has  1/R  to  ground  of  0 

C  to  ground  of  20 

1/L  to  ground  of  0 

From  node  2  to  node  3 

1/R  =  1 

C  =  1 

1/L  =  0 

From  node  2  to  node  4 

1/R  =  0 

C  =  0 

1/L  =  0 

From  node  2  to  node  5 

1/R  =  0 

C  =  0 

1/L  =  0 
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From  node 


2  to  node  6 
1/R  =  0 

C  =  0 

1/L  =  0 

Node  3  has  1/R  to  ground  of  0 

C  to  ground  of  250 
1/L  to  ground  of  0 

From  node  3  to  node  4 
1/R  =  0 

C  =  0 

1/L  =  78700 

From  node  3  to  node  5 
1/R  =  0 

C  =  0 

1/L  =  0 

From  node  3  to  node  6 

1/R  =  0 

C  =  0 

1/L  =  0 

Node  4  has  1/R  to  ground  of  0 

C  to  ground  of  0 

1/L  to  ground  of  0 

From  node  4  to  node  5 

1/R  —  2000 

C  =  0 

1/L  *  0 

From  node  4  to  node  6 

1/R  =  0 

C  =0 
1/L  =  0 

Node  5  has  1/R  to  ground  of  0 

C  to  ground  of  250 

1/L  to  ground  of  0 

From  node  5  to  node  6 

1/R  =  1 

C  =  1 

1/L  -  0 

Node  6  has  1/R  to  ground  of  20 

C  to  ground  of  20 

1/L  to  ground  of  0 

Source  admittance  =  1000000 

Lowest  frequency  considered  0. 100E-02 

No.  of  decades  considered  6 


APPENDIX  4 


This  appendix  presents  the  results  of  a  complete  response  analysis 
applied  to  the  same  circuit  as  in  Appendix  3. 

Complete  response  circuit  analysis  program,  version  1 

03-22-1985 

Complete  response  analysis  -for  a  circuit  with  6  nodes. 

time  step  =  .001  ,  number  of  steps  =  500 

node  1  has  1/R  to  ground  of  0 
C  to  ground  of  0 
1/L  to  ground  of  0 

From  node  1  to  node  2 
1/R  =  20 

C  =  0 

1/L  =  0 

From  node  1  to  node  3 
1/R  =  0 

C  =  0 

1/L  =  0 

From  node  1  to  node  4 
1/R  =  0 

C  =  0 

1  /L  =  0 

From  node  1  to  node  5 
1/R  =  0 

C  =  0 

1/L  =  0 

From  node  1  to  node  6 
1/R  =  0 

C  =  0 

1/L  =  0 

node  2  has  1/R  to  ground  of  0 

C  to  ground  of  20 
1/L  to  ground  of  0 

From  node  2  to  node  3 
1/R  =  1 

C  =  1 

1  /L  =  0 

From  node  2  to  node  4 
1/R  =  0 

C  =  0 

1/L  =  0 

From  node  2  to  node  5 
1/R  =  0 

C  =  0 


r,  r.  -r  s. 

■V-'.nV-  .V.s'V.V.'*  .wX- 


From  node 


2  to  node  6 
1/R  =  0 

C  =  0 

1/L  =  0 

node  3  has  1/R  to  ground  of  0 

C  to  ground  o-f  250 
1/L  to  ground  of  0 

From  node  3  to  node  4 

1/R  =  0 

C  =  0 

1/L  =  78700 

Initial  current  into  node  3  via  the  inductor  to  node  4  i 

From  node  3  to  node  5 

1/R  =  0 

C  =  0 

1  /L  =  0 

From  node  3  to  node  6 

1/R  =  0 

C  =  0 

1/L  =  0 

node  4  has  1/R  to  ground  of  0 
C  to  ground  of  0 
1/L  to  ground  of  0 

From  node  4  to  node  5 

1/R  -  2000 

C  =  0 

1  /L  =  0 

From  node  4  to  node  6 

1/R  =  0 

C  =  0 

1/L  =  0 

node  5  has  1/R  to  ground  of  0 

C  to  ground  of  250 
1/L  to  ground  of  0 

From  node  5  to  node  6 

1/R  =  1 

C  =  1 

1  /L  =  0 

node  6  has  1/R  to  ground  of  20 
C  to  ground  of  20 
1/L  to  ground  of  0 
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»‘*  •**  /■  «■*  •  •*"  «■-  •**  -  k  •  »  ■»  .  •  »  w  -  m 
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The  following  equivalencies  may  be  useful  in  applying  the  program.  In 
this  Appendix  we  use  impedances  (in  Ohms  and  Henries).  The  conversion  from  a 
Thevenin  to  a  Norton  Source  is. 


DISTRIBUTION  LIST 


Strategic  Defense  Initiative  Organization 
Directed  Energy  Weapons  Office 
The  Pentagon 

Office  of  the  Secretary  of  Defense 
Washington,  DC  20301-7100 

ATTN:  LTC  Richard  L.  Gullickson 

Commander 

Naval  Sea  Systems  Command 
Department  of  the  Navy 
Washington,  DC  20363 

ATTN:  NAVSEA/PMS  405  (Capt  R.L.  Topping) 
CDR  William  F.  Bassett 
Mr.  David  L.  Merritt 

Air  Force  Weapons  Laboratory  (NTYP) 

Kirtland  Air  Force  Base 
Albuquerque,  NM  87117 

ATTN:  LTC  James  Head 

U.S.  Army  Ballistics  Research  Laboratory 
Aberdeen  Proving  Ground,  MD  21005 

ATTN:  Dr.  D.  Eccleshall  (DRDAR-BLB) 

Lawrence  Livermore  Laboratory 
University  of  California 
Livermore,  CA  94550 

ATTN:  Dr.  R.  J.  Briggs 
Dr.  K.  Struve 
Dr.  W.  Barletta 
Dr.  D.  Prono 
Dr.  Y.P.  Chong 
Dr.  F.W.  Chambers 

Pulse  Sciences  Inc. 

14796  Wicks  Blvd. 

San  Leandro,  CA  94577 

ATTN:  Dr.  S.  Putnam 

Science  Applications  Inc. 

Security  Office 
5  Palo  Alto  Square,  Suite  200 
Palo  Alto,  CA  94304 

ATTN:  Dr.  R.R.  Johnston 
Dr.  Leon  Feinstein 


8.  Naval  Surface  Weapons  Center 
White  Oak  Laboratory 
Silver  Spring,  MD  20910 

ATTN:  Dr.  C.  M.  Huddleston,  R401 

Dr.  R.  B.  Fiorito,  R41 

Dr.  H.  S.  Uhm,  R41 

Dr.  Eugene  E.  Nolting,  H23 

9.  C.S.  Draper  Laboratories 
555  Technology  Square 
Cambridge,  MA  02139 

ATTN:  Mr.  E.  Olsson 

10.  Office  of  Naval  Research 
Department  of  the  Navy 
Arlington,  VA  22217 

ATTN:  Dr.  W.  J.  Condell  (Code  421) 

1 1 .  Avco  Everett  Research  Laboratory 
2385  Revere  3each  Pkwy. 

Everett,  MA  02149 

ATTN:  Dr.  Dennis  Reilly 

12.  Defense  Technical  Information  Center 
Cameron  Station 

5010  Duke  Street 

Alexandria,  Virginia  22314  (2  copies) 

13.  Naval  Research  Laboratory 
Washington,  D.  C.  20375 

ATTN:  T.  Coffey  -  Code  1001 
M.  Lampe  -  Code  4792 
M.  Friedman  -  Code  4700.1 
J.  R.  Greig  -  Code  4763  (50  copies) 

I.  M.  Vitkovit3ky  *  Code  4701 
W.  R.  Ellis  -  Code  4000 

S.  Ossakow,  Supt.  -  4700  (26  copies) 

Library  -  Code  2628  (20  copies) 

A.  Ali  -  Code  4700. IT 
D.  Book  -  Code  4040 

J.  Boris  -  Code  4040 

A.  Robson  -  Code  4760 
M.  Picone  -  Code  4040 
M.  Raleigh  -  Code  4763 
R.  Pechacek  -  Code  4763 
D.P.  Murphy  -  Code  4763 
R.F.  Fernsler  -  Code  4790 

J.  D.  Sethian  -  Code  4762 

K.  A.  Gerber  -  Code  4762 
G.  Joyce  -  Code  4790 

D.  Colombant  -  Code  4790 

B.  Hui  -  Code  4790 

R.  Hubbard  -  Code  4790 
Y.Y.  Lau  -  Code  4790 
Code  1220  -  1  copy 


Director  of  Research 

U.S.  Naval  Academy 

Annapolis,  MD  21402  (2  copies) 


15.  Mission  Research  Corp. 

1720  Randolph  Road,  S.E. 

Albuquerque,  NM  87106 

ATTN:  Dr.  Brendan  Godfrey 
Dr.  J.R.  Clifford 
Dr.  R.  Adler 
Dr.  G.  Kiuttu 

16.  McDonnell  Douglas  Research  Laboratories 
Dept.  223,  Bldg.  33,  Level  45 

Box  516 

St.  Louis,  MO  63166 

ATTN:  Dr.  J.C.  Leader 

Dr.  Evan  A.  Rose 

17.  Cornell  University 
Ithaca,  NY  14853 

ATTN:  Prof.  David  Hammer 

18.  Sandia  National  Laboratories 
Albuquerque,  NM  87185 

ATTN:  Dr.  Bruce  Miller,  1270 
Dr.  Carl  Ekdahl 
Dr.  M.  Mazarakis 
Dr.  C.  Frost 

19.  AFOSR/NP 

Bolling  Air  Force  Base,  Bldg.  410 
Washington,  DC  20331 

ATTN;  Capt.  H.  Pugh 

20.  University  of  Michigan 
Department  of  Nuclear  Engineering 
Ann  Arbor,  MI  48109 

ATTN:  Prof.  Terry  Kammash 

Prof.  Ronald  M.  Gilganbach 

21 .  SRI  International 
333  Ravenswood  Avenue 
Menlo  Park,  CA  94025 

ATTN:  Dr.  D.  Eckstrom 

22.  Los  Alamos  National  Laboratory 
Los  Alamos,  NM  87545 

ATTN:  Dr.  R.  Carlson 

Dr.  S.  Czuchlewski 


14.  Defense  Advanced  Research  Projects  Agency 
Directed  Energy  Office 
1400  Wilson  Blvd. 

Arlington,  VA  22209 

ATTN:  Dr.  Shen  Shey 


END 

p 

i 

i. 

FILMED 

■7i  i  >  ^  ■■ 

1-86 

DTIC 


